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Evcnt-by-event fluctuations and correlations between par- 
ticles produced in relativistic nuclear collisions are studied. 
The fluctuations in positive, negative, total and net charge are 
closely related through correlations. In the event of a phase 
transitions to a quark-gluon plasma, fluctuations in total and 
net charge can be enhanced and reduced respectively which, 
however, is very sensitive to the acceptance and centrality. If 
the colliding system experiences strong density fluctuations 
due, e.g., to droplet formation in a first-order phase tran- 
sition, all fluctuations can be enhanced substantially. The 
importance of fluctuations and correlations is exemplified by 
event-by-event measurement of the multiplicities of J/ty's and 
charged particles since these observables should anti-correlate 
in the presence of co-mover or anomalous absorption. 
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Event-by-event fluctuations have been measured at the 
SPS H§] and, with the higher multiplicities of RHIC 
and LHC, will become an important tool for studying 
the anomalous fluctuations and correlations that might 
remain following the phase transition to a quark-gluon 
plasma (QGP). Studies of event-by-event fluctuations at 
SPS energies do not indicate the presence of new physics 
H-||. It has, however, been proposed that large multi- 
plicity fluctuations can arise from density fluctuations or 
droplets created by a first-order phase transition Q . Re- 
cently, it has also been suggested that fluctuations in net 
charge should be suppressed in a QGP f|,@l§. Here, we 
shall consider multiplicity fluctuations in some general- 
ity to see how they are affected by conservation of total 
charge and strangeness, to understand the correlations 
between various measured fluctuations, and to show how 
their measurement can reveal interesting details of the 
collision process. 

Multiplicity fluctuations between various kinds of 
particles can be strongly correlated. As an example, con- 
sider the multiplicities of positive and negative pions, 
2V+ and A_, in a rapidity interval Ay for any relativis- 
tic heavy-ion experiment. (Similar analyses can be per- 
formed for any two kinds of distinguishable particles.) 
We define the fluctuation in any multiplicity A as 
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Empirically, the fluctuations of Eq. (|l|) are typically of 
order unity in relativistic nuclear collisions which is con- 
sistent with the expectations of Poisson statistics. 

The net positive charge from the protons in the collid- 
ing nuclei is much smaller than the total charge produced 
in an ultrarelativistic heavy-ion collision. For example, 
(N+) exceeds (A_) by only ~ 15% at in Pb+Pb collisions 
at SPS energies. The fluctuations in the number of posi- 
tive and negative (or neutral) pions are also very similar, 
ujn + — wat_ ■ Charged particle fluctuations have been es- 
timated in thermal as well as participant nucleon models 
H including effects of resonances, acceptance, and im- 
pact parameter fluctuations. By varying the acceptance 
and centrality, the degree of thermalization can actually 
be determined empirically || . Detailed analysis indicates 
that the fluctuations in central Pb+Pb collisions at the 
SPS are thermal whereas peripheral collisions are a su- 
perposition of pp fluctuations |l(| . 

The fluctuations in the total (N c h = N + + A_) and 
net (Q = N + — A_) charge are defined as 



((A+ ± A_) 2 ) - (N+ ± A_) 5 



(A+ 



•JV_) 
(N-) 



Ll> N _ ± C . 



(N ch ) 

where the correlation is given by 
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Fluctuations in positive, negative, total and net charge 
can be combined to yield both the intrinsic fluctuations 
in the numbers of N± and the correlations in their pro- 
duction as well as a consistency check. These quantities 
can change as a consequence of thermalization and a pos- 
sible phase transition. 

In practice, lon + ~ wjv_, so that the fluctuation in 
total charge simplifies to 



_ (N 2 ch ) - (N ch )' 



{Nch) 

and that for the net charge becomes 



uj n+ + C . 
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The fluctuation in net charge is related to the fluc- 
tuation in the ratio of positive to negative particles, 
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luq ~ (N + /N_)(N c h)uiN/N + /4: plus volume (or impact 
parameter) fluctuations |5||| . The virtue of this expres- 
sion is that volume fluctuations can in principle be ex- 
tracted empirically. Alternatively one can vary the cen- 
trality bin size or the acceptance H . In the following we 
shall assume that such "trivial" volume fluctuations have 
been removed. 

The analysis has so far been general and Eqs. (|[||) 
apply to any kind of distinguishable particles, e.g. pos- 
itive and negative particles, pions, kaons, baryons, etc. 
- irrespective of what phase the system may be in, or 
whether it is thermal or not. In the following, we shall 
consider thermal equilibrium, which seems to apply to 
central collisions between relativistic nuclei, in order to 
reveal possible effects on fluctuations of the presence of 
a quark-gluon plasma. 

Bosons/fermions have thermal fluctuations, lom = 
1 ± (np)/(n p ) where n p = (exp(e p /T) =F is the 

boson/fermion distribution function, which are slightly 
larger/smaller than those of Poisson statistics for a Boltz- 
mann distribution. Massless bosons, e.g. gluons, have 
uj b = £(2)/£(3) = 1.37 and massless fermions, e.g. 
quarks, have top = 2£(2)/3£(3) ~ 0.91 independent of 
temperature. Massive bosons have smaller fluctuations 
with, for example, w v = 1.11 and lj p = 1.01 when 
T = m n . 

In a thermal hadron gas (HG) as created in rela- 
tivistic in nuclear collisions, pions can be produced ei- 
ther directly or through the decay of heavier resonances, 
p, w, .... The resulting fluctuation in the measured num- 
ber of pions is 
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CON, — LON — — Itt^tt + fp^p + fuj^uj + 
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where f r is the fraction of measured pions produced from 
the decay of resonance r, and J2 r /r = 1< These mecha- 
nisms are assumed to be independent, which is valid in a 
thermal system. 

The heavier resonances such as p , u>, ... decay into 
pairs of tt + tt~ and thus lead to a correlation 
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Resonances reduce the fluctuations in net charge in a HG 
to ujq = 0.70 [|||]. In addition, total charge neutrality 
reduces fluctuations in net charge when the acceptance is 
large and thus increases correlations as will be discussed 
below. 

A phase transition to the QGP can alter both fluc- 
tuations and correlations in the production of charged pi- 
ons. To the extent that these effects are not eliminated by 
subsequent thermalization of the HG, they may remain 
as observable remnants of the QGP phase. As shown in 
Refs. ^,0], net charge fluctuations in a plasma of u, d 
quarks and gluons are reduced partly due to the intrinsi- 
cally smaller quark charge and partly due to correlations 
from gluons 



N f 
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where Nf is the number of quark flavors, g/ their charges, 
and N q the number of quarks. The total number of 
charged particles (but not the net charge) can be altered 
by the ultimate hadronization of the QGP. This effect 
can be estimated by equating the entropy of all pions to 
the entropy of the quarks and gluons. Since 2/3 of all 
pions are charged and since the entropy per fermion is 
7/6 times the entropy per boson in a QGP 



(N ch )c^((N g )+ 7 -(N q )), 
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where the number of gluons is (N g ) = (16/9Nf)(N q ). In- 
serting this result in (||) , we see that the resulting fluctua- 
tions are ujq — 0.18 in a two-flavor QGP (and ujq = 0.12 
for three flavors). As pointed out in 0, lattice results 
give loq ~ 0.25. However, according to || a substantial 
fraction of the pions are decay products from the HG, 
and the entropy of the HG exceed that of a pion gas by a 
factor 1.75 — 1.8. As described in the net charge fluc- 
tuations should be increased by this factor in the QGP, 
i.e. ujq ~ 0.33 in a two-flavor QGP, whereas it remains 
similar in the HG, loq ~ 0.6. 

However, if the high density phase is initially domi- 
nated by gluons with quarks produced only at a later 
stage of the expansion by gluon fusion, the production of 
positively and negatively charged quarks will be strongly 
correlated on sufficiently small rapidity scales. If, for ex- 
ample, the entropy density increases by an order of mag- 
nitude in going from a HG to QGP without additional 
net charge production, fluctuations in net charge will be 
reduced significantly, 
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The resulting fluctuation in net charge is necessarily 
smaller than that from thermal quark production as 
given by Eq. (@). A similar phenomenon occurs in string 
models where particle production by string breaking and 
qq pair production results in flavor and charge correla- 
tions on a small rapidity scale |ll]|l^] ■ The recently mea- 
sured charged particle density at midrapidity in central 
nuclear collisions at RHIC jl3| is only ~ 30 — 40% larger 
than pp scaled up by the nuclear mass as was also found 
at SPS energies jp. If entropy and multiplicities are pro- 
portional, the net and total charged particle fluctuations 
should be the same as at SPS according to Eq. (|Io|) un- 
less anomalous non-thermal fluctuations occur as will be 
discussed below. 

The strangeness fluctuation in kaons K might seem 
less interesting at first sight since strangeness is not sup- 
pressed in the QGP: The strangeness per kaon is unity, 
and the total number of kaons is equal to the number 
of strange quarks. However, if strange quarks are pro- 
duced at a late stage in the expansion of a fluid initially 
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dominated by gluons, the net strangeness will again be 
greatly reduced on sufficiently small rapidity scale. Con- 
sequently, fluctuations in net/total strangeness would be 
reduced / enhanced. 

The baryon number fluctuations have been estimated 
in a thermal model 0]. It is, however, not known how the 
annihilation of baryons and antibaryons in the hadronic 
phase affect these results. If only charged particles are 
detected, but not K°, K°, neutrons and antineutrons, 
the fluctuations have smaller correlations as compared 
to the total and net strangeness or baryon number. 

Total charge conservation is important when the 
acceptance Ay is a non-negligible fraction of the total 
rapidity. It reduces the fluctuations in the net charge as 
calculated within the canonical ensemble, Eqs. If 
the total positive charge (which is exactly equal to the 
total negative charge plus the incoming nuclear charges) 
is independently distributed according to the single par- 
ticle distributions, the resulting fluctuations within the 
acceptance are 



— 1 — facc , 
^Nch = 1 — face + 2 face 



(11) 
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where f acc = N^ ot J Ay (dN c h/dy)dy is the acceptance 
fraction or the probability that a charged particle falls 
into the acceptance Ay. Since charged particle rapidity 
distributions are peaked near midrapidity, charge con- 
servation effectively kills fluctuations in the net charge 
even when Ay is substantially smaller than the labora- 
tory rapidity, yi ao ~ 6 (11) at SPS (RHIC) energies. To- 
tal charge conservation also has the effect of increasing 
uj c h towards 2wjv + according to Eqs. ^ and jl^). Similar 
effects can be seen in photon fluctuations when photons 
are produced in pairs through p° — > 2j. In the WA98 
experiment, w 7 ~ 2 is found after acceptance corrections 
and eliminating volume fluctuations 

When the acceptance Ay is too small, particles from a 
thermal ensemble can diffuse in and out of the acceptance 
during hadronization and freezeout JtJ. Furthermore, 
correlations due to resonance production will disappear 
when the average separation in rapidity between decay 
products exceeds the acceptance. Each of these effects 
tend towards Poisson statistics when Ay^8y, where Sy 
denotes the rapidity interval that particles diffuse during 
hadronization, freezeout and decay. If luq is the canonical 
thermal fluctuation of Eqs. (0,^|), the resulting fluctuation 
after correcting for both Sy and total charge conservation 
is approximately 



exp 



Ay 



Ay + 2Sy 



26y 



Ay + 25y 



(1 - face) • (13) 



Here, the factor (1 — / acc ) is due to total charge conser- 
vation as in ( pi] ) . The remainder is fluctuations from two 
sources: a thermal one with fluctuations ojq and a ran- 
dom one with Poisson fluctuations, each weighted with 
the fraction of the charged particles they contribute. 



The resulting fluctuations in total and net charge are 
shown in Fig. 1 assuming u>n+ = ^ir — 1-1 and Sy — 
0.5. As mentioned above, f acc and Ay are related by 
the measured charge particle rapidity distributions [[J. 
Preliminary NA49 data plS agree well with the net and 
total fluctuations in a HG (C = 0.4) from Eqs. ( [l3|) and 
(M) . Residual volume fluctuations are significant for lon cH 
Q and have been estimated and subtracted. The curves 
apply to RHIC energies as well after scaling Sy with Ay. 
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FIG. 1. Acceptance dependence of thermal fluctuations in 
net charge (uoq of Eq. (|l3|), lower curves) and total (uiN ah , 
upper curves). Correlations increase from a hadron gas 
(C ~ 0.4), to a QGP (C ~ 0.8) and a pure gluon plasma 
(C ~ 1.1) (see text). 

The net charge fluctuations in a thermal HG corrected 
for finite acceptance and diffusion arc slightly below the 
value without any intrinsic correlations given by Eq. 
( |TT| ) . This reduction is due to the correlations in the HG 
that leads to ljq < 1 . In high energy pp collisions there 
are stronger rapidity correlations between unlike than 
like charged particles pT| , p^ | leading to similar magnitude 
for loq []. Therefore the net charge fluctuations does not 
vary by much going from peripheral pp-\ike high energy 
nuclear collisions to central collisions that are more likely 
to produce a thermal hadronic gas. The fluctuations in 
total charge are, however, very different because the total 
charge fluctuations in pp collisions increase dramatically 
with collision energy, uj^ h J c ~ 6 and ^YF^ — 20 as 
compared to uiff PS ~ 2.0 pd] , ^2[ . Peripheral collisions 
will therefore be very different from central ones and the 



1 The correlations are related to the two-body density dis- 
tributions, e.g., (N+N-) = J p^l_ (j/i, y2)dy\dy2, where the 
integral extends over 3/1,3/2 £ Ay. 
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centrality dependence should be studied carefully to as- 
sess the degree of thermalization before anomalous fluc- 
tuations due to phase transitions can be determined [||. 

Large non-thermal fluctuations can arise as a con- 
sequence of density fluctuations due, e.g. to droplet for- 
mation in first-order phase transitions. These could lead 
to large fluctuations in multiplicities [f||]]l4|] of charged 
particles and therefore also in the total and net charge. 
The above estimates for the fluctuations were of order 
unity. They implicitly assumed a uniform expanding 
system. Consider a scenario where the total multiplic- 
ity within the acceptance arises from a normal hadronic 
background component (Nhg) and from a second com- 
ponent (Nqgp) that has undergone a transition: 



N = N HG + N Q gp . 



(14) 



Its average is (N) = (Nhg) + (Nqgp)- Assuming that 
the multiplicity of each of these components is statisti- 
cally independent, the multiplicity fluctuation becomes 



— Uhg + (^QGP — UHG 



(Nqgp) 
(N) 



(15) 



Here, ujhg is the standard fluctuation in hadronic mat- 
ter u>hG — 1- The fluctuations due to the component 
that had experienced a phase transition, ujqgp, depend 
on the type and order of the transition, the speed with 
which the collision zone goes through the transition, the 
degree of equilibrium, the subsequent hadronization pro- 
cess, the number of rescatterings between hadronization 
and freezeout, etc. If thermal and chemical equilibration 
eliminate all signs of the transition, ujqgp — <-^hg- At 
the other extreme, the droplet scenario could produce 
wqgp ~ (N) ~ 10 2 — 10 3 if most hadrons arrive from a 
droplet so that either all or none fall into the acceptance 
Q. This is a promising signal worth looking for. Since 
droplets or density fluctuations are expected to be charge 
neutral, net charge fluctuations should vanish ujq ~ 
whereas u> c h — 2&n + ~ 2^qgp- 

General correlators between all particle species 
should be measured event-by-event, e.g., the ratios B 



(m/iNj) 
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u Nj (NiNj) - (Ni){Nj) 
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where iVjj are the multiplities in acceptances i and j of 
any particle type. (We assume that Nj is so large that it 
never vanishes.) In the presence of droplets, Ni and Nj 
would be strongly correlated in nearby rapidity intervals 
and at all azimuthal angles. 

As another example, consider correlations between 
multiplicities of charmonium particles N^ 7 ip = 
J/tyjip', x, ■■ and charged particle multiplicity (N) in a 
given rapidity interval Ay. If a ip is absorbed on co- 
movers or anomalously suppressed by QGP, one would 
expect anti- correlations because the number of co-movers 
and QGP should scale with the multiplicity N. By con- 
trast, direct Glauber absorption should not depend on 



the multiplicity of particles in Ay for a given central- 
ity since it is the result of collisions with participating 
nucleons in Glauber trajectories along the beamlinc. 

To quantify this anti-correlation, we model the absorp- 
tion of ips by simple Glauber theory 



i _ 



-lN/(N) 
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where N^ is the average number of ip's for given charge 
particle multiplicity N, and N® is the number of ips 
before co-mover or anomalous absorption sets in but 
after direct Glauber absorption on participant nucle- 
ons. In Glauber theory, the exponent is the prod- 
uct of the absorption cross section (o~ C ip), the absorber 
density (p c ), and the average path length (I) traversed 
in matter. The density, therefore also the exponent, 
is proportional to the multiplicity N with coefficient 
7 = — dlogN^/dlog N. In a simple co-mover absorption 
model for a system with longitudinal Bjorken scaling, 7 
can be calculated to be approximately Ea| 



E 



dN c (v C ipo-cip) 
dy AttR 2 



log(fl/r ) , 



(18) 



where dN c /dy is the co-mover rapidity density, a C ip the 
absorption cross section, v ci p the relative velocity, R the 
transverse size of the overlap zone, and To the formation 
time. Co- Mover absorption reduces the number of ip by a 
factor e -7 where 7 increases with centrality up to 7 ~ 1 
for typical parameters employed in co-mover absorption 
models @|. 

Since fluctuations in the exponent are small, 
j^/ujn I ' (N) 1, the anti-correlation is 



(NN.^)-(N)(N^) 
TO 



(19) 



It is negative and proportional to the amount of co-mover 
and anomalous absorption. It vanishes when the absorp- 
tion is independent of multiplicity (7 = 0). The rapidity 
interval should not be less than the typical relative rapidi- 
ties between co-movers and the or the rapidity range 
of a droplet. 

Even in central heavy- ion collisions ip's are rarely pro- 
duced and so iV^ = 1 or 0. In both cases one should mea- 
sure dNch/dy and average over events with and without 
a ip separately. If comovers absorb the ip we expect that 
(dN c h/dy) is slightly smaller for the events with a ip than 
without, leading to the negative correlation in Eq. (|l9|). 

A quantitative assessment of the average suppression 
of ips due to co-mover absorption versus direct ?/;-nucleon 
absorption has been debated ever since the first mea- 
surements of J/W suppression. The anticorrelations of 
Eq. (|l^) directly quantify the amount of co-mover or 
anomalous absorption and can therefore be exploited to 
distinguish between these and direct Glauber absorption 
mechanisms. In that respect it is similar to the elliptic 
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flow parameter for ip [jl5). For a sample of Nf„ rJ „ t the 
statistical uncertainty in 7 as determined by Eq. (y_9|) is 

~ \J (N) /^events ■ ^ we t a ke a rapidity bin Ay ~ 1 and 
consider central heavy ion collisions, (N) will range from 
~ 10 2 - 10 3 in going from SPS to RHIC energies. A 
sample of 10 4 — 10 5 i/js would be sufficient to determine 
7 with an accuracy of ±0.1. The analysis of the kaon to 
pion ratio by NA49 obtains such an accuracy by compar- 
ing to a mixed event analysis which removes systematic 
errors 0. 

The impact parameter fluctuations and correlations 
may not cancel exactly for the tp/N c h ratio, as they do 
for the 7r + /ir~ ratio jq], because their production mech- 
anisms differ. It is hard for the ip and soft for most of 
the charged particles and thus their multiplicities scale 
approximately with the number of binary NN collisions 
and the number of participants respectively. The number 
of NN collisions increase more rapidly with centrality and 
nuclear mass number (cx A 4 / 3 . As a result the impact pa- 
rameter correlations between the ip and N in Eq. ( Jigj ) 
will be slightly larger than the impact parameter fluctu- 
ations implicit in the second term in Eq. (|l6|). The dif- 
ference is a finite fraction of the total impact parameter 
fluctuations and will be of similar magnitude but oppo- 
site sign as correlations from co-mover absorption. The 
net impact parameter fluctuations will, however, depend 
on centrality and decrease as bin-size of the centrality 
cut decreases which should make it possible to separate 
it from other correlations. It would reveal independent 
information on the soft vs. hard production mechanisms. 

In summary, we have given a detailed analysis of to- 
tal and net charge fluctuations and correlations includ- 
ing total charge conservation and diffusion effects and 
how they depend on the acceptance. The correlations 
may increase if a QGP is formed resulting in a reduc- 
tion/enhancement of net/total charge by up to an order 
of magnitude depending on the model. It is important to 
measure fluctuations and correlations for various accep- 
tances as well as versus centrality and/or beam energy. 
At SPS energies the fluctuations in net charge actually 
increase slightly with centrality JlC| ] due to thermaliza- 
tion which is opposite to the predicted decrease in net 
charge fluctuations if a QGP is formed in central heavy- 
ion collisions. At RHIC and LHC energies the correla- 
tions vs. centrality will be much larger because the total 
charge fluctuations in pp collisions are w^ HIC ~ 6 and 
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uj^c ~ 20 as compared to u}f£* ~ 2.0 |11|,|12[. Periph- 
eral collisions will therefore be very different from cen- 
tral ones and the centrality dependence should be stud- 
ied carefully to assess the degree of thermalization before 
anomalous fluctuations due to phase transitions can be 
determined Q . It is important to measure all multiplicity 
fluctuations and correlations and understand the physical 
effects discussed above before a possible small increase in 
correlations can be attributed to the formation of a QGP. 
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